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The effects of mass-varying neutrinos on cosmic microwave background (CMB) anisotropies and 
large scale structures (LSS) are studied. In these models, dark energy and neutrinos are coupled 
such that the neutrino masses are functions of the scalar field playing the role of dark energy. We 
begin by describing the cosmological background evolution of such a system. It is pointed out 
that, similar to models with a dark matter/dark energy interaction, the apparent equation of state 
measured with SNIa can be smaller than -1. We then discuss the effect of mass-varying neutrinos on 
the CMB anisotropies and the matter power spectrum. A suppression of power in the CMB power 
spectrum at large angular scales is usually observed. We give an explanation for this behaviour and 
discuss different couplings and quintessence potentials to show the generality of the results obtained. 
We perform a likelihood analysis using wide-ranging SNIa, CMB and LSS observations to assess 
whether such theories are viable. Treating the neutrino mass as a free parameter we find that the 
constraints on the coupling are weak, since CMB and LSS surveys give only upper bounds on the 
neutrino mass. However, fixing a priori the neutrino masses, we find that there is some evidence 
that the existence of such a coupling is actually preferred by current cosmological data over the 
standard ACDM cosmology. 
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I. INTRODUCTION 

Recent cosmological observations indicate that the ex¬ 
pansion of the universe is accelerating 0 - 0 - It follows 
from General Relativity that the dominant energy com¬ 
ponent today must have negative pressure. Many candi¬ 
dates have been proposed over the years, including scalar 
field models, which are well motivated from the point of 
view of particle physics theories, see e.g. 0-(l]. The 
main prediction of these types of models is that the dark 
energy equation of state becomes a dynamical quantity, 
and can vary from the usual value of w = — 1 for a cosmo¬ 
logical constant. Although such models are very attrac¬ 
tive, they are plagued with several theoretical difficulties, 
such as the stability of the potential under quantum cor¬ 
rections Q or why the dark energy scalar field seems not 
to mediate a force between normal matter particles [sj. 
In addition, the energy scale of the scalar held is put in 
by hand and usually not connected to a more fundamen¬ 
tal energy scale. However, attempts have been made to 
address these problems, such as models with ultralight 
pse udo-Nambu-Goldstone bosons (see, for example, }10| . 
{ill and [0; for a review, see 0). 

It is expected that any explanation for dark energy will 
involve physics beyond the standard model of particle 
physics. Recently, a new class of models have been pro¬ 
posed, which entertain the idea of a possible connection 


between neutrinos and dark energy. Their theoretical 
and observational consequences have already been stud¬ 
ied very extensively |l4|-[48|. The main motivation for a 
connection between dark energy and neutrinos is that the 
energy scale of dark energy (C>(10~ 3 ) eV) is of the order 
of the neutrino mass scale. In these models the neutrino 
mass scale and the dark energy scale are linked to each 
other, and hence the observed non-zero neutrino masses 
(see |0-[5lj) cannot be understood without an under¬ 
standing of dark energy. Also, one may hope that these 
models might provide an explanation for the coincidence 
problem [181 ]. 

In this paper we investigate the cosmology of neutrino 
models of dark energy. We take into account the full evo¬ 
lution of the neutrinos, i.e. studying the relativistic and 
non-relativistic regimes and the transition in between. 

Armed with a complete numerical model for the evo¬ 
lution of the coupled neutrinos, we compare the back¬ 
ground evolution with Supernova data and study how 
the modified perturbations affect the cosmic microwave 
background radiation (CMB) temperature anisotropies 
and large scale structures (LSS) matter power spectrum. 
We thereby present the details of the results outlined in 
0 and discuss other forms of coupling between dark 
energy and neutrinos and potentials for the dark energy 
field.' 

The paper is organized as follows: In Section II 
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we discuss the background evolution of the coupled 
dark energy-neutrino system in the context of a typical 
quintessential potential. In Section III we derive the evo¬ 
lution equations for cosmological perturbations in neu¬ 
trino models of dark energy, and present the modified 
CMB and matter power spectra. In Section IV we discuss 
other couplings and potentials, such as inverse power-law 
potentials and field-dependent couplings. In Section V 
we compare our theory with data, using a public Markov- 
Chain Monte-Carlo data analysis program. We conclude 
in Section VI. 


II. THE COSMOLOGICAL BACKGROUND 
EVOLUTION 

In a flat, homogeneous, Friedmann-Robertson-Walker 
universe with line-element 

ds 2 = a 2 (r) (—dr 2 + 5ijdx l dx 4) , (1) 

the Einstein equations describe the evolution of the scale 
factor a(r): 

*-(§)’-*<*» < 2 > 

s(£) = < 3 > 

In these equations, p(r) and p(r) are the total energy 
density and pressure respectively and the dot refers to 
the derivative with respect to conformal time r. Defin¬ 
ing fli = pi/p c , where p c is the critical energy density for 
a flat universe and pi are the energy densities of the in¬ 
dividual matter species, the equations above require that 
O = - Oj = 1. In the following we will set 87 tG = 1. 

In our model we consider a universe with the usual 
energy-matter composition. At early times, the energy 
density is dominated by the relativistic species - radia¬ 
tion and highly relativistic neutrinos. As the universe ex¬ 
pands the energy density in radiation decays, and the uni¬ 
verse becomes matter dominated. The dominant matter 
species is assumed to be Cold Dark Matter (CDM), which 
is non-relativistic, weakly interacting and behaves like a 
perfect pressureless fluid. At this time there are also 
contributions to f 1 from baryons and neutrinos (which 
having cooled behave in a manner similar to CDM). At 
late times the matter energy densities also decay away, 
and we enter the dark energy dominated epoch. In com¬ 
mon with standard quintessence models we describe the 
dark energy sector using a dynamical scalar held with 
potential U (</>), where the form of the potential is chosen 
(and fine-tuned) to produce the necessary late time ac¬ 
celeration. The energy density and pressure of the scalar 
held are dehned by the usual expressions, 

P<t> = 2^ 2 + ( 4 ) 

= ^2</> 2 - Vffl. (5) 


In this paper, we consider the consequences of a cou¬ 
pling between neutrinos and dark energy. To describe 
this coupling, we follow 0: the coupling of dark energy 
to the neutrinos results in the neutrino mass becoming a 
function of the scalar held, i.e. m v = and so the 

mass of the neutrinos changes as the scalar held evolves. 
For our purposes it does not matter if the neutrinos are 
Majorana or Dirac particles, and for simplicity we as¬ 
sume three species of neutrinos with degenerate mass 1 It 
is well known {§] that the light mass of the quintessence 
potential results in it being highly unstable to radiative 
corrections, and that the addition of a coupling between 
the dark energy and other matter species only serves to 
further exacerbate this problem. In this regard it is im¬ 
portant that both the quintessence potential and the neu¬ 
trino mass are regarded as classical, effective quantities, 
which already include radiative corrections. 

It is also important to note that our theory differs sig¬ 
nificantly in one key aspect from the work of [18]. In 
our models, the dark energy sector is described by a light 
scalar field, with a mass which is at most of order H. 
The potential chosen by Fardon et al. was such that the 
mass of the scalar field is much larger than H for most 
of its history, and this can have significant implications 
upon the behaviour of the neutrino background and the 
growth of perturbations [H] as we will discuss later. 

To fully describe the evolution of cosmological neu¬ 
trinos, we must calculate their distribution function 
/ (x l ,p l , r) in phase space. An important fact to note 
is that even though the neutrinos interact with dark en¬ 
ergy, we treat the interaction classically and, as will be 
shown in eq. (P2121) . they can be thought as free-falling in 
a metric given by 

9a0 = m v {4>fg a p. ( 6 ) 

Thus, the theory we are going to consider is a special 
type of scalar-tensor theory, in which the scalar degree 
of freedom couples only to neutrinos. It follows that the 
neutrino phase-space density is incompressible and we 
can treat the neutrinos as collisionless particles through¬ 
out the period of interest as long as we keep track of 
the evolution of the neutrino mass. We therefore need to 
solve the Boltzmann equation in collisionless form simul¬ 
taneously with the scalar held evolution equations. Once 
the distribution function is known, the pressure and en¬ 
ergy density of the neutrinos can be calculated. In this 
Section we will discuss the background evolution only; in 
the next Section we will discuss cosmological perturba¬ 
tions in these models. 


1 In fact, such an assumption is quite natural and has no strong 

consequences to this work. In the mass regions detectable in 

astronomical observations, the three neutrino masses are nearly 
degenerate. Adding to that, cosmology is in leading order sensi¬ 

tive to y~l mi/. 
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The energy density stored in the neutrinos is given by 


Pi> = “4 J q 2 dqdnef 0 (q), (7) 

and the pressure by 

Pv = 3 a 1 / ( 8 ) 


where fo(q) is the usual unperturbed background neu¬ 
trino Fermi-Dirac distribution function 


/o(e) 


g s 1 

h 3 P e e / kBT ° + 1 ’ 


(9) 


and e 2 = q 2 + m 2 (4>)a 2 (<7 denotes the comoving mo¬ 
mentum). As usual, g s , lip and kp stand for the num¬ 
ber of spin degrees of freedom, Planck’s constant and 
Boltzmann’s constant respectively. In the following we 
will assume that the neutrinos decouple whilst they are 
still relativistic, and therefore the phase-space density 
only depends upon the comoving momentum. Taking 
the time-derivative of eq. 0, it can be easily shown that 


n 1T") r rri 

p v + 3 H {p„ + p v ) = —— -if {p„ - 3 p v ) . ( 10 ) 

dq> 

We describe the dark energy sector using a scalar field 
with potential V(<f >). Taking into account the energy con¬ 
servation of the coupled neutrino-dark energy system, 
one can immediately find that the evolution of the scalar 
field is described by a modified Klein-Gordon equation 

cj) + 2Hcf + a — = -a 2 ——— {p v - 3p „). (11) 

dcp dip 


This equation contains an extra source term with respect 
to the uncoupled case, which accounts for the energy ex¬ 
change between the neutrinos and the scalar field. 

For the remainder of this Section and the next, we 
consider a typical exponential form for the dark energy 
potential, namely 


V(<t>) = V 0 e~^ ( 12 ) 


and define er 


A We also choose to take 


m„{4>) = M o e 0(t> 


(13) 


The exponential potential can produce a non-scaling 
solution that may give late time acceleration, depen ding 
upon the steepness of the potential, er (see e.g. [H, [35^ 
HH, 113: HI) • In an uncoupled system with er < \/6 there 
exists a critical point that is stable for er 2 < 3(1 + w), 
where w stands for the equation of state of matter or 
radiation, and in which f = 1. This solution will lead 
to acceleration provided that er < y/2. The existence of 
scaling solutions depends upon the equation of state of 
the other components present in the universe. Choosing 


er 2 > 3(1 + w) leads to a scaling solution with = 
3(1 + w)/a 2 0]. (See also HU, who use the exponential 
potential as a dark energy model.) The requirement that 
the present day dark energy density is f~ 0.7 is hard 
to reconcile with the scaling solution at early times, since 
in this case it follows that f2^ = 4/er 2 , whilst big bang 
nucleosynthesis requires that the dark energy density in 
the early universe is very small [38| . 

In this Section we focus our discussion on models with 
a < v 2 j which with an appropriate choice of Vo, can 
provide late-time acceleration with f2^ ~ 0.7 today (note 
that the late-time attractor f 2 ^ = 1 lies in the future). 
This choice of a also ensures that the energy density in 
the form of dark energy at the time of BBN is very small, 
because for early times the quintessence field is frozen and 
acts like a cosmological constant with an energy density 
similar to the observed dark energy density today. 

The presence of neutrino coupling can potentially af¬ 
fect this result, as the coupled held begins to evolve at 
earlier redshifts (z ~ 10 '), however the fraction of the 
total energy stored in the scalar held at early times re¬ 
mains insignificant. Note that this choice of potential 
reduces to the cosmological constant case for a perfectly 
hat potential with zero coupling. 

From the neutrino energy conservation equation m, 
and for our choice of m v {4>) and one can see that 

the dynamics of the scalar held can be described by an 
effective potential 

V eS = V{cf>) + (p v - 3 p v ) e?*, (14) 

where p v = p„e - ^ and p„ = p v e~ 0 ^ are independent 
of <f>. It can be shown that the effective potential will 
only have a minimum when /3a > 0. For the neutrinos 
we have numerically evaluated the integrals 0 and 0 , 
which then have been used in the Klein-Gordon eq. m 
to hnd the evolution of the scalar held. 

Figure 1 shows some typical examples of how the cou¬ 
pling of the neutrinos to the scalar held causes the mass 
of the neutrinos to evolve with time. Deep within the 
radiation dominated epoch, at times when the neutrinos 
are highly relativistic, the scalar held is Hubble damped 
and therefore the neutrino mass is (almost) constant. For 
quintessence models, <f> is at most of order H , and there¬ 
fore for relativistic species the coupling term in equation 
(fTUl) is clearly suppressed relative to the Hubble damp¬ 
ing term. As the universe expands the neutrinos cool and 
become non-relativistic at a temperature corresponding 
to the neutrino mass. Hence, the extra coupling terms 
in equations 01 and m become more and more im¬ 
portant, allowing energy to be exchanged between the 
neutrinos and the scalar held. This interaction causes 
the scalar held, and hence the neutrino mass, to evolve. 

It is important to note the two different types of be¬ 
haviour seen in Figure 1 for the evolution of the neu¬ 
trino mass. For models which have (3a > 0 the effective 
potential possesses a minimum, and after some time the 
held passes through this minimum, slows down, stops and 
eventually rolls back towards the minimum. For models 
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FIG. 1: Plot showing the evolution of the neutrino mass for an 
exponential potential and coupling (solid line: (3 = 0, A = 1; 
short dashed line: /3 = 1, A = 1; dotted line: (3 = —0.79, 
A = 1; long dashed line (3 = 1, A = O.5.). In all models, we 
have arranged that m v = 0.314 eV today. 

which do not posses an effective minimum (f> is always 
negative, and the scalar field will continue to roll down 
the effective potential unimpeded. We compare the be¬ 
haviour of our light scalar field with the heavy acceleron 
field used in fl8j] - in their model the scalar field sits in 
the effective minimum of its potential for most of the 
time during the cosmic history, and it is the evolution of 
the effective minimum which drives the dynamics of the 
neutrino mass. As discussed in [l8|, the mass of their 
neutrinos increases as the universe expands, whereas in 
our model the neutrinos are heavier in the past and be¬ 
come lighter (although as we will discuss in Section IV, 
suitable choices of coupling and potential can realize cou¬ 
pled neutrino-quintessence models with neutrinos which 
are lighter in the past). 

For the model described in this section the coupled 
neutrinos are heavier in the past than uncoupled neu¬ 
trinos, which implies that the energy density stored in 
the neutrinos is larger than would normally be expected. 
This means that the evolution of the neutrino density 
parameter depends on the evolution of the neutrino 
mass, which in turn depends on the choice of the cou¬ 
pling (3 and the slope of the potential a. This can be 
seen from Figure 2. The coupling of neutrinos to dark 
energy significantly alters the evolution of the cosmolog¬ 
ical background. In particular it can be seen that the 
extra energy stored in the neutrinos in the past can alter 
the redshift of matter radiation equality. 

Typically one expects non-relativistic neutrinos to be¬ 
have in a similar manner to CDM, however the interac¬ 
tion between the neutrinos and the scalar field modifies 
the scaling behaviour of the non-relativistic neutrinos. 
It can be seen that the neutrino energy density dilutes 
away faster than that of CDM, which is especially notable 
for large values of coupling. The evolution of <f> caused 




?+/ 
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FIG. 2: Background evolution: In the upper panel, we plot 
the evolution of the density parameters for a model with f3 = 
0, A = 1. In the middle panel the corresponding plot with (3 = 
1 is shown, whilst the lower plot shows the case f3 = —0.79, 
A = 1. (Neutrinos: solid line, CDM: dot-dashed line, scalar 
field: dotted line and radiation: dashed line.) In all cases, the 
mass of the neutrinos is = 0.314 eV today. We consider 
a flat universe with = 0.022, 0l c h 2 = 0.12, f2„/i 2 = 0.01 
and h = 0.7. 


by the transfer of energy between the coupled neutrinos 
and scalar field also results in the energy density of the 
quintessence field becoming dominated by kinetic energy. 
Finally at a redshift of the order of unity the potential en¬ 
ergy of the scalar field begins to dominate, and all other 
matter species decay away. 

A final point we would like to raise is the fact the 
apparent equation of state measured by an observer is 
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not given by 



(15) 


where p$ and p$ are defined in eqns (0D and (0). One of 
the usual assumptions made in the measurements of the 
dark energy equation of state using supernovae is that 
matter (dark, baryonic or neutrinos) is decoupled from 
dark energy. At low redshifts, all these components are 
assumed to scale like a~ 3 . This is clearly not the case 
with the coupled neutrinos here. It was pointed out in 
[a 3 |. that the apparent equation of state is given by 2 


w ap 


W 4 , 

1 — x ’ 


(16) 


with 

Pv, 0 


In this equation, the subscript 0 denotes the quantities 
at the present epoch. We emphasise that this quantity is 
not the effective equation of state of dark energy, which 
is defined as 


m v ((t) 0 ) 


- 1 


(17) 


p<fr + 317/90 (1 + Wett) = 0, (18) 


while assuming that the neutrino density and neutrino 
pressure evolve according to Eq. m- Using the Klein- 
Gordon equation, one finds that the effective equation of 
state can be written as 


/ 3(j) p„ 

w eS = w </> + — —. 

3 H p^ 


(19) 


In Fig. 3, we plot the apparent equation of state u> ap 
as a function of redshift. Note that w ap can be less than 
— 1, as was pointed out in [ssj . [56j and [54| in the context 
of models with dark matter/dark energy interaction. 

To conclude this part, even if dark energy couples only 
to a subdominant component such as neutrinos, the ap¬ 
parent equation of state can be less than —1, without 
introducing phantom fields. As it can be seen from Fig¬ 
ure 3, the apparent equation of state might even cross 
the boundary w = —1. 


III. PERTURBATION EVOLUTION 



FIG. 3: The apparent equation of state, as defined in eq. 
m as a function of redshift 3 (solid line: /? = 0, A = 1; 
short dashed line: /3 = 1, A = 1; dotted line: (3 = —0.79, 
A = 1; long dashed line f3 = 1, A = 0.5.) 


(Fo r a review of cosmological perturbation theory, see 
S, @ or [HO 

As already mentioned in the last section, to fully de¬ 
scribe the evolution of cosmological neutrinos, we must 
calculate their distribution function / (a:*, p*,t) in phase 
space. We can treat the neutrinos as collisionless parti¬ 
cles throughout the period of interest, and hence we can 
find the neutrino distribution function by solving the col- 
lisionless Boltzmann equation |57] 


d£,d^dl,dqd£ Add/ ( , 

dr dr dx i dr dq dr dn i ’ 

where the comoving momentum is qt = api. It is con¬ 
venient to rewrite the comoving momentum in terms of 
its magnitude and direction: q l = qri 1 . The last term 
in equation 1HD is a second order quantity and will be 
neglected in the following linear perturbation formalism. 

The path of a neutrino in spacetime is governed by the 
general action 


S = - j m iy (4 > )\/~ds 2 1 


( 22 ) 


Let us now turn our attention to the evolution of cos¬ 
mological perturbations in our model. We adopt the con¬ 
ventions of Ma and Bertschinger [57( and work in the 
synchronous gauge, taking the line element to be 

ds 2 = —a 2 dr 2 + a 2 (Sij + hij) dx l dx^. (20) 


which can be minimised to derive the neutrino geodesic 
equation 3 


P° 


dpp 

!h 


T p p a pP 

OifJ 


, d In m„ dip 
dip dx p 


(23) 


where is the proper momentum of the neutrino. Tak¬ 
ing the zeroth component of this equation and using the 


2 The authors of [H3l called this quantity an effective equation 
of state. However, we will define an effective equation of state 

below. 3 Please note Errata at the end of this paper. 
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relation P° = ea 2 , one finds that the comoving three- 
momentum of the neutrinos is given by 3 

^ qhijmrij . (24) 

This equation does not depend explicitly on the coupling 
or the scalar field perturbations. Following [57], we write 
the phase space distribution of the neutrinos as a zeroth 
order distribution plus a small perturbation 

f(x\pj,r) = f 0 (q) [l + ® ( x\q,n,j,T )] . (25) 


Substituting this expression into the Boltzmann equation 
and performing a Fourier transformation, we find 3 


d^f 

77 


i— (k • n) T 

e 

d In /q [. 
ding V 


h + 6f] 
2 


(k-n) 2 


= 0 


(26) 


In this equation (and in eqs. (l29l) and (l3ll) . given later), rj 
and h are the standard scalar parts of the metric pertur¬ 
bation hij. It is clear that equation (PZ(I1) does not contain 
terms proportional to dlnm„/<I</>. Therefore, the equa¬ 
tions for the neutrino hierarchy derived in [571 ] do not 
change 3 . However, the expressions for the perturbed neu¬ 
trino energy density and neutrino pressure, which will be 
calculated using /, are modified. The perturbed energy 
density is given by 


r 1 f 2 , , m f t dlnm„ m 2 a 2 \ 
dp v = J q dqdflfo ( e'k + ——-— 5<j>J (27) 

which can be written as 

S P V = ~4 J g 2 dg ^£/o(g)^+^ C?1 ^™' y {pv-Zpv)- (28) 


Similarly, the expression for the perturbation in the neu¬ 
trino pressure is given by 


Sp v = 


1 

3a4 


q 2 dq dilfoiq) 


q 2 T „ , d\n.m v q 2 m 2 cd 

— % ~ S( P—T7 - 5 

e dcp e J 


The expressions for the neutrino shear and energy flux 
remain unchanged as they do not depend explicitly upon 
m v . Finally, the perturbed Klein-Gordon equation is 
given by: 


5<p + 2H5(j> - 

din m u 

d(j) 


,2 2 d 2 V 

1 + »v 

(dp,, - 3 8p v ) 


Scj) + —hcj) = 


(29) 


d 2 In to j, 
d(j) 2 


5cj)(p v - 3 p v ) 


To calculate the temperature anisotropy spectrum and 
matter power spectrum we apply these modifications 
to CAMB [61]. This code calculates the linear cosmic 
background anisotropy spectra by solving the Boltzmann 



10 100 1000 
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FIG. 4: The CMB anisotropy spectrum (unnormalized) for 
exponential coupling and potential. Solid line: (3 = 0, A = 1; 
short-dashed line: (3 = 1, A = 1; dotted line: (3 = —0.79, 
A = 1; long-dashed line: f3 = 1, A = 0.5. Error bars denote 
WMAP data. 
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FIG. 5: Evolution of the sum of the metric perturbations 
<f> + Solid line: (3 = 0, A = 1; short-dashed line: (3 = 1, 

A = 1; dotted line: (3 = —0.70, A = 1; long-dashed line: 
(3 = 1, A = 0.5. The scale is k = 10 -3 Mpc -1 . 


equation which governs the evolution of the density per¬ 
turbations, and integrating the sources along the pho¬ 
ton past light cone. To ensure the accuracy of our cal¬ 
culations, we directly integrate the neutrino distribu¬ 
tion function, rather than using the standard velocity 
weighted series approximation scheme. We do not con¬ 
sider lensing effects, nor tensor contributions. 

The results of the neutrino-dark energy coupling on 
the temperature anisotropy spectrum can be seen in Fig¬ 
ure 4. The most obvious modifications to the anisotropy 
spectrum occur for large angular scales, with ( < 100, al¬ 
though for some choices of parameters the positions and 
relative heights of the peaks are also affected. We gener- 
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FIG. 6: Plot of the matter power spectrum. From the top 
curve to the bottom curve: (/3 = 0, A = 1), (/3 = 1, A = 0.5), 
( j3 = —0.79, A = 1). The matter power spectrum for (/3 = 1, 
A = 1) is indistinguishable from the (/3 = 0, A = 1) curve. 


ally observe an increase in power for 10 < £ < 100 , whilst 
for £ < 10 we find either an excess or reduction in power 
depending upon our choice of parameters. Note that this 
is in marked contrast to models of coupled CDM, where 
an increase in power on large scales is usually observed 
[63 |. For the models where the neutrinos were much heav¬ 
ier in the past than today, we also observe a slight shift 
in the acoustic peaks and a change in their relative am¬ 
plitudes. 

For scales larger than a degree (£ < 100), the dominant 
contribution to the anisotropy spectrum is the Integrated 
Sachs-Wolfe Effect (IWS). This arises due to the evolu¬ 
tion of the gravitational potentials along the photon path 
from the surface of last scattering. The modification to 
the cosmological background arising from the neutrino 
coupling can have a significant effect upon the evolution 
of the perturbations. In particular there is a larger energy 
density in neutrinos in coupled models during the transi¬ 
tion period when the neutrinos become non-relativistic. 
As a result, the intermediate regime between radiation 
and matter domination is prolonged, and so the evolution 
of the gravitational potentials are significantly modified. 
The evolution of the sum of the metric perturbations $ 
and U/ is shown in Figure 5. The modifications to the 
behaviour of the metric perturbations for the different 
models is immediately apparent. 

For very large scales (£ < 20) anisotropies arise pri¬ 
marily from the late time Integrated Sachs-Wolfe effect 
(ISW), which is caused by the evolution of the metric 
perturbations for redshifts in the range 0 < z < 2. In 
particular, and p„ as well as the equation of state of 
dark energy affect the late time behaviour of cosmological 
perturbations. As mentioned above, the evolution of the 
scalar field is influenced by the presence of the coupling 
to the neutrinos and hence the equation of state of dark 
energy depends upon (3. Likewise, the clustering prop¬ 


erties of dark energy depends on the neutrino coupling 
(see [63| for a discussion on the clustering of dark energy 
and its impact on the CMB). The neutrinos will gener¬ 
ally tend to fall into the potential wells of dark matter, 
although at a rate slightly dependent on the coupling to 
the scalar field. The scalar field itself will cluster together 
with the neutrinos and thereby affecting the gravitational 
potential. 

Let us turn our discussion to the evolution of neutrino 
perturbations. Figure 6 shows the effects of neutrino cou¬ 
pling on the matter power spectrum. Here we typically 
observe damping, and our results appear similar to stan¬ 
dard models of CDM and hot dark matter, where a sim¬ 
ilar reduction in power could be achieved with a heavier 
neutrino mass. 

We can use the perturbed part of the energy momen¬ 
tum conservation equation for the coupled neutrinos 


rpll 

7; M 


d In nriv 
d(f> 




(30) 


to calculate the evolution equations for the neutrino per¬ 
turbations (T a a stands for the trace of the neutrino en¬ 
ergy momentum tensor and the semicolon denotes the 
covariant derivative). Taking 7 = 0 we derive the equa¬ 
tion governing the evolution of the neutrino density con¬ 
trast , S u = whilst taking 7 = * (spatial index) yields 

the velocity perturbation equation 6 v = ikiV l v , with the 
coordinate velocity v l v = dx l /dr: 


S v = 3 H 


(h + P4>) ( 


w v - S v — (l+w„) ( 6 V + ^ 


d/3 


+ (3 (1 — 3w v ) 54> + -j7<t>S<t> (1 — 3 w v ) ■ 

dtp 


(31) 
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w v 1 + w v 
v k 2 84> — (3{\ — 3wv)4>0v — k 2 a„. (32) 


The variable oy represents the neutrino anisotropic stress 
and we have used the more general definition (3 = 
dlnmv/dcj), which in general might be not constant. 
Furthermore the neutrino equation of state is given by 
W„ = Pv/Pv 

It is the presence of the additional coupling terms in 
these expressions for the growth of the neutrino density 
and velocity perturbations, as well as the modifications to 
the evolution of the cosmological background, which al¬ 
ters the behaviour of the neutrino perturbations in com¬ 
parison with the standard uncoupled case. 

Figure 7 shows the evolution of the neutrino and CDM 
density contrasts, comparing the uncoupled model with 
an extreme coupled case with (3 = —0.79, for which the 
mass of the neutrinos is m v « 2.5 eV at 2 > 1400 but 
To;, = 0.3 eV today. Deep inside the radiation dominated 
epoch the neutrinos are highly relativistic and their den¬ 
sity contrast grows in a similar manner to radiation. In 
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FIG. 7: Evolution of the neutrino (lower curves) and CDM 
(upper curves) density contrasts. The solid line shows the 
uncoupled case, i.e. j3 = 0, whereas the dashed line shows the 
case for (3 = —0.79. The scale is k = O.IMpc -1 . 


the uncoupled models the growth of the density pertur¬ 
bations of the neutrinos makes a transition to matter -like 
behaviour once the neutrinos become non-relativistic. 
The neutrinos will fall into the CDM potential wells, 
which is the dominant component at recombination. At 
small wavelengths (large wavenumbers), the neutrinos 
undergo freestreaming, which prevents neutrinos from 
clustering at an arbitrary small scale. The freestreaming 
length scale after the neutrinos become non-relativistic 
can be estimated to be [H3] (reinstating Newton’s con¬ 
stant G) 


fcf s (a) 


47 t Gpa 2 

«med 


(33) 


where p is the background total density. The median 
neutrino speed is given by: 

u med = 15a _1 ^ km s_1 - (34) 


Since the neutrino momentum decays like a -1 , the 
neutrino velocity behaves like (am v (a ))~ 1 , taking into 
account that the neutrino mass evolves with time. 
Freestreaming stops as soon as k < kf s , allowing the 
neutrino density contrast to grow. This behaviour can 
clearly be seen in Figure 7 for both the uncoupled and 
coupled cases. At around 2 ss 10 4 the neutrinos become 
non-relativistic and start to freestream immediately, as 
can be seen from the oscillating behaviour of 5„. At this 
stage k{ s < k. However, as soon as k[ s = k freestream¬ 
ing stops, and S v can grow unimpeded. The case with 
neutrino-dark energy coupling differs from the uncoupled 
case since in the result shown the neutrinos are heavier in 
the past, so for a given redshift kf s is larger. This means 
that freestreaming stops earlier than in the uncoupled 
case. This behaviour is apparent in Figure 7 (dashed 


lines), where we see that S u starts to grow earlier than in 
the uncoupled case (solid lines). The neutrino-coupling 
also has an effect on the growth of the density contrast 
itself since we observe that S u grows slower than in the 
uncoupled case. This is probably because the rate of 
gravitational infall of the neutrinos tends to be reduced 
by the presence of the much less clustered dark energy. 
Also, the coupling has a slight effect on the growth of the 
dark matter density contrast, which arises from the fact 
that the background evolution is modified. 


IV. ANOTHER CHOICE OF COUPLING AND 
POTENTIAL 

So far we have restricted our discussion to one choice of 
quintessence-neutrino coupling and one form for the dark 
energy potential. At this stage, the reader might wonder 
whether the results obtained so far are simply due to 
our choice of potential and coupling. For a scalar field 
with standard kinetic term, the exponential potential is 
not a favored model for a quintessential potential, since 
the initial value of the scalar field has to be fine-tuned to 
obtain scalar field domination today [HU . The interesting 
alternative possibility of a global attractor unfortunately 
does not seem to be viable due to the large perturbation 
growth {62]. The coupling of the scalar field to neutrinos 
does not cure the fine-tuning problem of the exponential 
potential. 

For our second form of coupling we choose 

m v {(j)) = Moe^ 2 , (35) 

which was also recently used in a model with dark mat¬ 
ter/dark energy coupling in (65]. The effect of this choice 
is that the coupling function dh\m v /dfi becomes field- 
dependent, whereas it has been constant so far. Depend¬ 
ing on how the field evolves with time, the coupling can 
either grow or become smaller during the cosmic history. 
Field-dependent couplings are not uncommon in higher 
dimensional theories and can appear in brane-world the¬ 
ories, see for example [66] or [67l |. 

For the potential, we choose an inverse power-law 
potential, which is a well-motivated candidate for a 
quintessence field (see e.g. 168] and [69]). To be con¬ 
crete, we use 



With these choices for the potential and coupling, the 
effective minimum will exist if /3 > 0. 

The results for the neutrino-mass evolution, the ap¬ 
parent equation of state and the CMB anisotropy power- 
spectrum are shown in Figures 8, 9 and 10. 

The biggest difference to the case of a purely expo¬ 
nential coupling is that for positive /? the neutrinos are 
lighter in the past, as can be seen from Figure 8. Thus, 
with a convenient choice of potential and coupling, the 
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FIG. 8: The evolution of the neutrino mass for the inverse 
power-law potential with m„(<(>) = rrioe^ (solid line: (3 = 0; 
short dashed line: (3 = 0.2; dotted line: (3 = —0.2; long dashed 
line /? = 0.27.) 


FIG. 10: The CMB anisotropy spectrum (unnormalized) for 
an exponential potential and m v ((j)) = moe^ (solid line: 
(3 = 0; short dashed line: /3 = 0.2; dotted line: (3 = —0.2; 
long dashed line (3 = 0.27.) 



with negative /3, in which a reduction of power at low 
multipoles can be observed. The cases with positive (3 
can not be distinguished from the uncoupled case. The 
reason is that the neutrino density is smaller in the past 
than in the uncoupled case for this choice of potential 
and coupling. Hence, neutrinos are less important for 
the dynamics of the universe, and their imprint upon the 
CMB is correspondingly reduced. 

In essence, the physical explanations of the model pre¬ 
sented in our earlier paper [28| and in Sections II and III 
remain valid even for other choices of the potential and 
couplings, since they show how to relate the general be¬ 
haviour of a dynamical neutrino mass to the cosmological 
evolution. 


FIG. 9: The evolution of the apparent dark energy equation 
of state for the inverse power-law potential with m„(0) = 
moe^ (solid line: /3 = 0; short dashed line: /3 = 0.2; dotted 
line: /3 = —0.2; long dashed line (3 = 0.27.) 


neutrinos can become heavier as the universe expands. In 
the case of a negative (3, the effective potential does not 
have a minimum and the neutrinos become lighter as the 
universe expands. The results for the apparent equation 
of state are shown in Figure 9. The results are similar to 
the ones found in Section II: w ap can be smaller than — 1 
and can cross the boundary of the cosmological constant 
with w = —1. As it can be seen in Figure 9, the appar¬ 
ent equation of state varies substantially in the redshift 
range 2 = 0 — 2 if (3 is non-zero. A strongly, however, 
varying equation of state is not preferred by the data. 
Finally, the effects on the CMB anisotropies are similar 
to the ones found in Section III, as can be seen from Fig¬ 
ure 10. The only visible deviation from /3 = 0 is the case 


V. PARAMETER CONSTRAINTS 

In the earlier sections we demonstrated that models 
of coupled dark energy and neutrinos could produce a 
detectable signature in cosmological surveys. Indeed, the 
modifications to the background evolution (mainly to the 
dark energy equation of state), temperature anisotropy 
spectrum and matter power spectrum should allow us to 
constrain our model using current data sets. 

We perform our likelihood analysis using CosmoMC 
0 - This program uses a Markov-Chain Monte-Carlo 
(MCMC) engine to efficiently explore the cosmological 
parameter space. Typically we run five chains for each 
simulation, with no less than 35,000 samples per chain. 
We perform the usual convergence checks on the indi¬ 
vidual chains to ensure that the chains have fully sam¬ 
pled our parameter space. As well as visually confirming 
that the individual chains converge, we check the Gelman 
and Rubin R statistic (variance of chain means/mean of 
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chain variances) for each parameter, and ensure that the 
Raftery and Lewis convergence diagnostic is satisfied. 

We take advantage of the wide range of cosmologi¬ 
cal data which is currently available to constrain our 
model. The CMB temperature anisotrop y sp ectrum is 
constrained using WMAP [n| [z2| , CBI [z|, ACBAR 
[Z3 and VSA [75} datasets. The neutrino-dark energy 
coupling can also affect the formation of large scale struc¬ 
ture which is sensitive to the neutrino mass, and so we 
use data from the Sloan Digital Sky Survey [76} to fur¬ 
ther constrain our model. Data from the Supernova Cos¬ 
mology Project can also be used to constrain the 
equation of state of dark energy, and thus place further 
constraints on our model. 

We choose to perform the data analysis using our usual 
choice of potential and coupling, namely V{<j>) = Voe~ a ^ 
and m v ((j>) = We choose to focus on these poten¬ 

tials because they embody the typical behaviour observed 
for most models of coupled neutrinos, and they easily re¬ 
duce to the standard ACDM case ((3 = a = 0). These 
potentials also have the advantage that the initial choice 
of 4>i does not affect the evolution of the cosmology, as 
changes to the initial choice of <f> are equivalent to re¬ 
scalings of the mass parameters Mq and Vo- For general 
choices of potentials and couplings this useful degeneracy 
does not exist, as the neutrino-dark energy coupling can 
severely restrict the range of attractor solutions. Conse¬ 
quently the increased number of fine-tuned free param¬ 
eters required for these models would compromise the 
goodness of fit compared to simpler models. 

Throughout our analysis we assume a flat universe, 
with fi = 1. Initially we use the standard parameteriza¬ 
tion for our cosmological model, and vary the following 
parameters: Df,h 2 , SIcdm^ 2 , h , z re (the redshift of reioni¬ 
sation), fl„/i 2 , n s (the spectral index), 10 lo A s (the initial 
scalar perturbation amplitude) and the dark energy pa¬ 
rameters a and (3. We show the results from this initial 
analysis is Figure [IT] 

Clearly Figure [TT] shows that current cosmological data 
places no constraints on our new coupling parameters. It 
is well known that the current best fit analysis of cosmo¬ 
logical data can only place an upper limit on the mass 
of the neutrino, and that massless neutrinos are not ex¬ 
cluded by most cosmological data sets. However the re¬ 
sults presented so far in this paper require that the neu¬ 
trinos having a significant mass; indeed for low mass neu¬ 
trinos the effects observed in this paper become largely 
insignificant. Despite this, as we will show later, the 
strength of the neutrino coupling can still be constrained 
by the requirement for later time acceleration. 

Current cosmological data requires that the universe 
contains approximately 70% of dark energy, 30% of dark 
matter and some minor quantities of baryons and neu¬ 
trinos. In the dynamics of our model there could exist 
only one critical point able to guarantee such propor¬ 
tions in its vicinity (see e.g. jb2 ) • When exactly reached, 
this point is characterized by the total domination of the 
scalar field and exists only for |oj < y/6 « 2.5. Indeed 


FIG. 11: Posterior constraints for 9 parameter model de¬ 
scribed in the text. Shading denotes the mean likelihood of 
the samples, whilst the contours show the 68% and 95% con¬ 
fidence limits from the marginalised distribution. Solid lines 
on the 1-D plots show the marginalised posterior, whilst the 
dotted curves denote the mean likelihood of the samples. 

our computed likelihood will be shown to be contained 
in such boundaries. A preferred range for the parame¬ 
ter /3 is more difficult to predict, however, because of the 
previously described effects induced on the background 
and perturbation at different stages of the evolution, al¬ 
though it is clear that very large values of (3 will not be 
favoured by the data. But it is very important to empha¬ 
size that values of (3 of the order of unity are perfectly 
acceptable. 

Another consideration when choosing the scalar field 
parameters is the BBN constraint resulting from the early 
time modification of the scale factor evolution due to the 
presence of dark energy [H, (38). As discussed earlier, the 
neutrino coupling to dark energy in our model does not 
modify the energy density of the dark energy at the time 
of BBN (when the neutrinos are highly relativistic, and 
the coupling terms are negligible). For a quintessence 
potential with a scalar field dominated late time attrac¬ 
tor the BBN constraints of fH} are easily satisfied. We 
have confirmed this numerically for our coupled mod¬ 
els (for instance, in the case of (3 = 1, 12 ^ ~ 10 -25 at 
z ~ 10 8 ). Note that, in principle, our parameter space in¬ 
cludes both the late-time attractor and scaling solutions 
(that evolve like radiation or matter and which would 
provide a non-negligible contribution to the energy den¬ 
sity at the time of BBN). We find, however, that the 
scaling solution for the exponential potential is already 
strongly disfavored by the observations that we have used 
for our analysis and therefore additional BBN constraints 
would not modify our findings. 

Recent works [78} have used data from X-ray clusters 
to reduce the uncertainties on the lower bound for the 
neutrino mass, finding a value for the neutrino masses 
of Yl 171 ” = 0-56tQ;26 e ^- Notice though that the up¬ 
per bounds could in fact go up to X) m v = 2 eV, de¬ 
pending on the datasets used and the assumed priors 
|79(. Measurements of atmospheric neutrino oscillations 
suggest that there is at least one neutrino species with 
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Parameter 

Mean 

likelihood 

68% 

interval 

95% 

interval 

Ocdm/i 2 

0.102 ±0.004 

0.099 - 0.106 

0.096 - 

0.110 

Zre 

17.6 ±3.7 

16.0 - 19.8 

10.7- 

23.2 

a 

0.43 ±0.32 

0.29 - 0.60 

0.13 - 

0.95 

p 

0.75 ±0.64 

0.64 - 0.98 

0.11 - 

1.18 

n s 

0.96 ±0.01 

0.95 - 0.97 

0.93 - 

0.99 


TABLE I: Marginalised parameter constraints for our 6 pa¬ 
rameter model with fixed m „ = 0.2 eV, Q,bh 2 = 0.022 and 
h = 0.72. For this model we find x 2 /dof = 1570.1/1459. This 
compares with a x 2 /dof = 1610.1/1461 for a best-fit ACDM 
model using the same parameter set with a = /3 = 0. 


Parameter 

Mean 

likelihood 

68% 

interval 

95% 

interval 

Dcdm ^! 2 

0.100 ±0.003 

0.097 - 0.104 

0.094 - 

0.107 

Zre 

20.5 ±3.1 

19.4 - 22.0 

15.0 - 

25.1 

a 

0.52 ±0.29 

0.40 - 0.67 

0.00 - 

0.97 

P 

0.62 ±0.21 

0.58 - 0.74 

0.15- 

0.86 

n s 

0.97 ±0.01 

0.96 - 0.99 

0.95 - 

1.00 


TABLE II: Marginalised parameter constraints for our 6 pa¬ 
rameter model with fixed m „ = 0.3 eV, Q,bh 2 = 0.022 and 
h. = 0.72. In this case we find x 2 /dof = 1593.7/1459. This 
compares with a x 2 /dof = 1636.8/1461 for a best-fit ACDM 
model using the same parameter set with a = /3 = 0. 


m v > 0.05 eV [49]. The Heidelberg-Moscow experiment, 
searching for neutrino-less double beta decay claim detec¬ 
tion of electron neutrino mass m„ e between 0.2 eV and 
0.6 eV, with best fit m„ e = 0.36 eV [8Gj. 

We therefore choose to perform our analysis using two 
values of the neutrino mass today, m„ = 0.2 eV and 
TOj, = 0.3 eV, to investigate whether models of neutrino- 
dark energy coupling could in principle be constrained 
if neutrinos were independently confirmed to have a sig¬ 
nificant (m u > 0.1 eV) mass, consistent with current 
experiments measuring the neutrino mass. 

By choosing to fix the value of the neutrino mass today, 
we are required to specify the current value for the en¬ 
ergy density stored in neutrinos and value of the Hubble 
constant as the neutrino mass, the critical energy density 
in neutrinos and the Hubble constant are related via the 
usual formula 


f Ivh 2 


93.2 eV' 


(37) 


We also choose to fix the value of £l b h 2 as we do not ex¬ 
pect the behaviour and constraints of the baryon energy 
density to be significantly modified by our neutrino-dark 
energy coupling, as the observed effects on the anisotropy 
spectrum are largely limited to relatively low multipoles. 

The values for Hq and Q b h 2 can be determined inde¬ 
pendently from the cosmological data used in our MCMC 
analysis. The value of Hq = 72 km s _1 Mpc -1 can be 
obtained from the best fit of the Hubble Space Telescope 
Key Project [81], whilst the baryon density parameter 
fl b h 2 = 0.022 is favoured by Big Bang Nucleosynthesis 
models [Hj. We are therefore left with a cosmological 
model requiring 6 parameters: /3, er, Hcdm^ 2 , z re , A s 
and n s . 

The results of the MCMC analysis for neutrinos with 
a mass today of m u = 0.2 eV can be found in Table QJ 
whilst the results for rn v = 0.3 eV are given in Table 
El We quote the marginalised probability distributions 
and confidence intervals. Figures fl2l fe [T3l show the 2D 
probability distributions for the m„ = 0.2 eV and m u = 
0.3 eV models respectively. 

As expected, the neutrino coupling has little effect on 
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FIG. 12: Posterior constraints for 6 parameter model, with 
m v = 0.2 eV, Q. b h 2 = 0.022 and h = 0.72. Shading denotes 
the mean likelihood of the samples, whilst the contours show 
the 68% and 95% confidence limits from the marginalised dis¬ 
tribution. Solid lines on the 1-D plots show the marginalised 
posterior, whilst the dotted curves denote the mean likelihood 
of the samples. 


the value of com h 2 as the peak structure of the tem¬ 
perature anisotropy spectrum is largely unaffected. The 
value found for n s is also within the usual range for pa¬ 
rameter analysis of cosmological models which neglect 
tensor contributions. 

For both the m v = 0.2 eV and m v = 0.3 eV models we 
find that non-zero values of neutrino coupling strengths 
are preferred by the data. We also see that for these 
models a non-zero value for a is preferred over the usual 
cosmological constant, although a = 0 is not excluded at 
the 68% level. This is not surprising since we have seen 
that the equation of state of the dark energy for some 
choices of parameters in our coupled models is entirely 
consistent with the preferred value of w ap ~ — 1 found 
from supernova surveys. It is clear that models with 
heavier neutrinos allow stronger constraints to be placed 
upon the strength of the coupling. Indeed, for the 0.3 eV 
neutrinos we find that neutrino-dark energy coupling is 
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FIG. 13: Posterior constraints for 6 parameter model, with 
m v = 0.3 eV, Qbh 2 = 0.022 and h = 0.72. Shading denotes 
the mean likelihood of the samples, whilst the contours show 
the 68% and 95% confidence limits from the marginalised dis¬ 
tribution. Solid lines on the 1-D plots show the marginalised 
posterior, whilst the dotted curves denote the mean likelihood 
of the samples. 


preferred at the 1 sigma confidence level. This is to be 
expected as a larger neutrino mass today is equivalent to 
a higher energy density in neutrinos, and so any modifi¬ 
cation to the neutrino evolution will have a larger impact 
on the CMB and LSS for models with greater densities 
of neutrinos. In particular we have seen that there exist 
a range of non-zero (3 values capable of reducing power 
at low CMB multipoles. This last effect is most probably 
the cause for the relative peak in the likelihood for (3 of 
order unity. Furthermore a sharper drop at large values 
of (3 is observed in the likelihood most likely to limit the 
excessive growth at multipoles 10 < i < 100. 

It is however important to make clear that these con¬ 
straints rely upon the assumption that the neutrino mass 
is known, and that the neutrinos have a mass m v > 
0.1 eV. Although this assumption is consistent with cur¬ 
rent neutrino experiments, we can only make the state¬ 
ment that should the mass of the neutrino be found to 
be greater than 0.1 eV, then cosmological data can be 
used to constrain the strength of any neutrino-dark en¬ 
ergy coupling; indeed we find that there is some evidence 
that the existence of such a coupling is actually preferred 
by current cosmological data over the standard ACDM 
cosmology. 


VI. CONCLUSIONS 

We have investigated models of dark energy which cou¬ 
ple a quintessence scalar field to massive neutrinos. In 
these models, dark energy and neutrinos are coupled such 
that the neutrino masses become functions of the scalar 
field. The effects of such models on the cosmological 


background evolution, on the cosmic microwave back¬ 
ground anisotropies, and on the formation of large scale 
structures were analyzed. Additionally, we have also per¬ 
formed a likelihood analysis on the parameter space of 
such theories. 

We have focused on two specific models: In the first, 
the coupling between neutrinos and dark energy is con¬ 
stant and the quintessential potential is an exponential. 
The second model, which is better motivated from the 
particle physics point of view, has a neutrino-coupling 
which depends on the quintessence field (hence changes 
with time), whilst the scalar field has a power-law poten¬ 
tial. In spite of some specific differences between these 
two models (such as the energy density stored in the 
scalar field at recombination for example), the effects of 
the coupling on the CMB anisotropies and on the matter 
power spectrum are nevertheless explainable by the ba¬ 
sic mechanisms that we have identified earlier. Namely, 
the coupling modifies the background history and induces 
an ISW contribution to the CMB spectrum; the matter 
power spectrum is modified by the magnitude of the neu¬ 
trino mass during structure formation. Given the gener¬ 
ality of these explanations, the conclusions drawn from 
this investigation could probably be applied to any sim¬ 
ilar model with a neutrino-dark energy coupling. 

It is important to note that in our models, the dark 
energy sector is described by a light scalar field, with a 
mass which is at most of order H. This is in contrast 
to previous models [l|j] in which the mass of the scalar 
field is much larger than H for most of its history. The 
latter can have significant effects upon the behaviour of 
the neutrinos and the growth of their perturbations, and 
which is difficult to reconcile with current astronomical 
data [HI. 

Solving the collisionless Boltzmann equation for the 
neutrinos, we have investigated the relativistic and non- 
relativistic regimes and the transition period in between. 
Initially the neutrinos are highly relativistic, and dur¬ 
ing this period the quintessence field is frozen. The 
mass of the neutrinos therefore remains constant. As 
the neutrinos become non-relativistic they begin to ex¬ 
change energy with the quintessence field via the coupling 
term. At a temperature scale comparable to the neu¬ 
trino mass, the neutrinos become non-relativistic, whilst 
the quintessence field is dominated by kinetic energy. It 
is at this point that the neutrino mass begins to evolve 
significantly. The details of this behaviour and evolution 
depends on the choice of the coupling (3 and the potential 
parameter a. In fact, the masses of the neutrinos can be 
heavier or lighter in the past depending on the choice of 
potential and coupling parameters. 

The coupling of neutrinos to dark energy slightly al¬ 
ters the evolution of the cosmological background. It was 
found that similarly to models with a dark matter/dark 
energy interaction, the apparent equation of state mea¬ 
sured with Type la Supernovae at high redshift can be 
smaller than —1, without introducing phantom fields, 
and might even cross the boundary w = — 1. 
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The most obvious modifications to the CMB 
anisotropy spectrum occur for large angular scales, with 
t < 100, where the dominant contribution to the 

anisotropies is generated by the Integrated Sachs-Wolfe 
Effect (IWS). This arises due to the evolution of the 
gravitational potentials along the photon path from the 
surface of last scattering. The modification to the cos¬ 
mological background arising from the neutrino coupling 
can also have a significant effect upon the evolution of 
the perturbations. We generally observe an increase in 
power for 10 < t < 100, whilst for i < 10 we find ei¬ 
ther an excess or reduction in power depending upon our 
choice of parameters. For the models where the neutri¬ 
nos were much heavier in the past than today, we also 
observe a slight shift in the peaks and a modification in 
their relative amplitude. 

The matter power spectrum exhibits free-streaming 
damping even in the presence of dark energy-neutrino 
coupling. However, since the damping scale is mainly 
dependent on the value of the neutrino mass at the end 
of their relativistic stage, our results appear similar to 
the standard models with CDM and hot dark matter in 
which the mass is fixed at the relativistic plateau. It is 
obvious that the mass infered from the damping of the 
matter power spectrum is, in general, different from the 
neutrino mass measured with experiments in the labora¬ 
tory. 

We performed a likelihood analysis using SNIa, CMB 
and LSS datasets. Initially, we used the standard pa¬ 
rameterization for our cosmological model, characterized 
by exponential dependence of the dark energy potential 
and neutrino mass on the scalar field. For a flat uni¬ 
verse we varied all of the matter parameters, the Hubble 
constant, the initial power spectrum spectral index and 
amplitude and the instantaneous reionization parameter 
z le . As expected, the cosmological data did not place 
strong constraints on our new coupling parameters. This 
is no surprise, since it is well known that the current best 
fit analysis of cosmological data can only place an upper 
limit on the mass of the neutrino, and a zero neutrino 
mass is not excluded by most cosmological data sets. An 
interesting outcome was that couplings of order unity are 
perfectly acceptable with the actual data. 

To proceed, we chose to perform an analysis using two 
values of the neutrino mass today, to„ = 0.2 eV and 
m u = 0.3 eV, to investigate whether models of neutrino- 
dark energy coupling could in principle be constrained if 
neutrinos were independently confirmed to have a signif¬ 
icant mass (rn v >0.1 eV), consistent with current exper¬ 
iments. 

For both the rn v = 0.2 eV and m v = 0.3 eV models we 
found that non-zero values of neutrino coupling strengths 
of order unity are preferred by the data. We also saw that 
for these models a non-zero value for a is preferred over 
the usual cosmological constant, although a = 0 is not 
excluded at the 68% level. Models with heavier neutrinos 
allow stronger constraints to be placed upon the strength 
of the coupling. Indeed, for the 0.3 eV neutrinos we found 


that neutrino-dark energy coupling is preferred at the 1 
sigma confidence level. 

One should note that these constraints rely upon the 
assumption that the neutrino mass is known, and that 
the neutrinos have a mass m v >0.1 eV. Although this 
assumption is consistent with current neutrino experi¬ 
ments, we can only make the statement that should the 
mass of the neutrino be found to be greater than 0.1 eV, 
then current cosmological data can be used to constrain 
the strength of any neutrino-dark energy coupling. 
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Erratum 

The original version of th is p aper contained a typo 
and a mistake, as noted by [83]. This section contains 
the corrected equations and Figures, as published in our 
Erratum [84 ]. 

The geodesic equation (23) contains a typo and should 
read 


We are grateful to K. Ichicki for correspondence. The 
correct equations (2)-(4) have been derived in [83). 


P [ 


,dPP 

dr 


VP' = 


d In m v d(j) 
d<f> dx p 


(38) 


A subtle error occurred in eq. (24), which should read 


dq 

dr 


--qhijmrij - a 2 3 4 


m 2 dlnrriv d<j> dx l 
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(39) 


It follows that a scalar field dependent term should be 
included in the Boltzmann equation (26), giving: 
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(40) 


Therefore, the dipole equation for the neutrino hierarchy 
derived in [57] is subject to a change represented by a 
new term once again dependent on the scalar held: 



a 


FIG. 14: Evolution of the metric variables <F + T. Solid line: 
/3 = 0, A = 1; short-dashed line: /3 = 1, A = 1; dotted line: 
P = —0.70, A = 1; long-dashed line: p = 1, A = 0.5. The 
scale is k = 10~ 3 Mpc -1 . 


1 q, 


<Fi = --fc(4- 0 - 2<F 2 ) - --k 

O 6 O 6 


1 q a 2 m 2 d In m v d In /o 


q 2 d(j> ding 


Scj) 
(41) 

This modification will have an effect on the ISW effect, 
which is less pronounced than reported in our paper. The 
corrected evolution of the metric variables $ + 4' is shown 
in Figure [HI 

This will effect the anisotropies in the CMB, which are 
shown in Figure fl5l 

On the other hand, we do not find changes to the mat¬ 
ter power spectrum or the evolution of the neutrino den¬ 
sity contrast at the scale given in Fig. 7 in our origi¬ 
nal paper. At smaller scales we register small differences 
in the neutrino and scalar held density contrasts, which 
leads to the mentioned changes in the ISW effect. Note 
that the background evolution reported in |28|, Us) is not 
affected. 

We do not attempt to redo the comparison of our 
model with data, as the new WMAP 3-year data have 
been published since [86}, and these models where anal¬ 
ysed in [83]. 



10 100 1000 

Multipole l 


FIG. 15: The CMB anisotropy spectrum (unnormalized) for 
exponential coupling and potential. Solid line: (3 = 0, A = 1; 
short-dashed line: p = 1, A = 1; dotted line: p = —0.79, 
A = 1; long-dashed line: P = 1, A = 0.5. 
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